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Abstract 

A measurement of charm and beauty dijet photoproduction cross sections at the ep collider 
HERA is presented. Events are selected with two or more jets of transverse momentum 

p\ 1<2) > 11(8) GeV in the central range of pseudo-rapidity —0.9 < r/ etl < 2 ) < 1.3. The 
fractions of events containing charm and beauty quarks are determined using a method 
based on the impact parameter, in the transverse plane, of tracks to the primary vertex, 
as measured by the HI central vertex detector. Differential dijet cross sections for charm 
and beauty, and their relative contributions to the flavour inclusive dijet photoproduction 
cross section, are measured as a function of the transverse momentum of the leading jet, 
the average pseudo-rapidity of the two jets and the observable x° bs . Taking into account 
the theoretical uncertainties, the charm cross sections are consistent with a QCD calcula- 
tion in next-to-leading order, while the predicted cross sections for beauty production are 
somewhat lower than the measurement. 
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1 Introduction 



A measurement is presented of charm and beauty production in ep collisions at HERA using 
events with two or more jets at high transverse momentum. The measurement is carried out 
in the photoproduction region in which a quasi-real photon, with virtuality Q 2 ~ 0, is emitted 
from the incoming positron and interacts with a parton from the proton. Differential charm and 
beauty dijet cross sections are measured and compared to calculations in perturbative quantum 
chromodynamics (pQCD) performed to next-to-leading order (NLO). 

In pQCD calculations, the photoproduction of charm and beauty proceeds dominantly via 
the direct photon-gluon fusion process jg — > cc or bb, where the photon interacts with a gluon 
from the proton to produce a pair of heavy quarks in the final state. Previous charm mea- 
surements have confirmed this prediction [1-4]. In leading order QCD models a successful 
description of the data is obtained when additional contributions from processes involving re- 
solved photons are taken into account [2, 3]. In such resolved photon processes the quasi-real 
photon fluctuates into a hadronic state before the hard interaction and thus acts as a source of 
partons. In the massless scheme, a large fraction of these resolved photon processes is due to 
heavy quark excitation, in which one of the partons that enters the hard interaction is a heavy 
quark (c or b) originating from the resolved photon or the proton. 

In this analysis events containing heavy quarks are distinguished from light quark events by 
the long lifetimes of c and b flavoured hadrons, which lead to displacements of tracks from the 
primary vertex. This technique, based on the precise spatial information from the HI silicon 
vertex detector, was introduced in recent HI measurements of the charm and beauty structure 
functions F% c and F 2 bb in deep inelastic scattering [5,6], and is now applied to dijet events in pho- 
toproduction. This analysis provides the first simultaneous measurement of charm and beauty in 
photoproduction, extending to larger values of transverse jet momentum than previous measure- 
ments [1-4,7-10]. The regions of small transverse momentum and large pseudo-rapidities are 
excluded from the measurement, due to trigger requirements and the limited angular acceptance 
of the vertex detector. 

The differential charm and beauty dijet cross sections are measured as functions of the trans- 
verse momentum of the leading jet pf 1 , of the mean pseudo-rapidity fj of the two jets, and of 
the variable a;° bs which, in a leading order picture, corresponds to the fraction of the photon's 
energy in the proton rest frame that enters the hard interaction. For direct photon-gluon fusion 
processes s° bs ~ 1, while for resolved photon processes x° bs can be small. The measured dif- 
ferential charm and beauty cross sections, together with the measured flavour inclusive cross 
sections, are used to determine the relative contribution from charm and beauty events to dijet 
photoproduction. The results are compared with calculations in perturbative quantum chro- 
modynamics at next-to-leading order and with predictions from Monte Carlo simulations in 
which leading order matrix elements are implemented, and contributions from higher orders are 
approximated using parton showers. 

This paper is structured as follows: In section the experimental apparatus is briefly de- 
scribed. Event and track selections are detailed in section The method to determine the 
contributions of charm and beauty events is outlined in sectional Theoretical calculations per- 
formed in the framework of perturbative QCD are discussed in section The cross section 
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measurements and their systematic uncertainties are presented in section Properties of a 
heavy quark enriched data sample are investigated in section |7] A summary of the results is 
given in section [5J 



2 HI Detector 

The HI detector is described in detail in [11]. Charged particles emerging from the ep in- 
teraction region are measured by the central tracking detector (CTD) in the pseudo-rapidity 
range —1.74 < rj < 1.74 1 . The CTD consists of two large cylindrical central jet drift cham- 
bers (CJCs), two z chambers and two multi-wire proportional chambers arranged concentrically 
around the beam-line in a magnetic field of 1.15 T. The CTD provides triggering information 
based on track segments from the CJC in the r-(p plane, transverse to the beam direction, and 
on the z position of the vertex from the multi-wire proportional chambers. The CJC tracks 
are linked with hits in the Central Silicon Tracking detector (CST) [12], which consists of two 
cylindrical layers of silicon strip sensors, surrounding the beam pipe at radii of 57.5 mm and 
97 mm from the beam axis. The detector provides hit resolutions of 12 /im in r-<p with an 
average efficiency of 97%. 

For CTD tracks with CST hits in both layers the impact parameter 5, i.e. the transverse 
distance of closest approach to the nominal vertex, can be measured with a resolution of ag ~ 
33 fim © 90 /im/p t [GeV]. The first term represents the intrinsic resolution and includes the 
uncertainty of the CST alignment, the second term corresponds to the contribution from multiple 
scattering in the beam pipe and the CST, which depends on the transverse momentum p t of the 
track. 

Charged and neutral particles are measured in the liquid argon (LAr) calorimeter which sur- 
rounds the tracking chambers and covers the range —1.5 < 77 < 3.4 and a lead-scintillating fibre 
calorimeter SpaCal, covering the backward region (—4.0 < i] < —1.4) [13]. The measurements 
from CTD and calorimeters are combined to reconstruct the final state particles [14]. The lu- 
minosity determination is based on the measurement of the Bethe-Heitler process (ep — > epj), 
where the photon is detected in a calorimeter located downstream of the interaction point in the 
positron beam direction. 



3 Event and Track Selection 

The data sample corresponds to an integrated luminosity of 56.8 pb -1 and was recorded with 
the HI experiment during the years 1999 and 2000. During this time HERA was operated 
with positrons of 27.5 GeV energy and protons of 920 GeV. The events were triggered by a 
combination of signals from the calorimeters, the central drift chambers and the multi-wire pro- 
portional chambers. Photoproduction events are selected by requiring that there be no isolated 
high energy electromagnetic cluster detected in the calorimeters consistent with a signal from 

'The pseudo-rapidity is given by r\ = — In tan(0/2), where 6 is measured with respect to the z-axis given by 
the proton beam direction. 
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the scattered positron. This restricts the range of negative four-momentum transfer squared to 
Q 2 < 1 GeV 2 . The inelasticity y is calculated using the hadronic final state [15], and the mea- 
surement is restricted to the range 0.15 < y < 0.8. The jets are reconstructed from the final 
state particles using the inclusive k t algorithm [16] in the p t recombination scheme [17], with 
distance parameter R — 1 in the -q-cj) plane. The event selection requires at least two jets in the 
central pseudo-rapidity range —0.9 < r\ < 1.3 with transverse energy p t 1(2) > 11(8) GeV. 

CTD tracks are selected which are linked to hits in both r-0 layers of the CST These tracks 
are required to have a transverse momentum above 500 MeV and a polar angle in the range 
30° < # trac k < 150°. Events are selected which contain at least one selected track associated to 
one of the two leading jets. The final sample consists of 80769 events. 

To correct for detector effects, such as detector resolutions and inefficiencies, large sam- 
ples of charm, beauty and light quark events are generated using the Monte Carlo program 
PYTHIA [18] (for details see section |5]). All samples are passed through a detailed simulation 
of the HI detector response based on the GEANT program [19] and the same reconstruction 
and analysis algorithms as used for the data. 

Figure[l] shows the distributions of the transverse momentum of the leading jet p t jetl , and 
of the jet with the second highest transverse momentum p^ 2 (figures [l^ and b), of the mean 
pseudo-rapidity fj = (rf etl + rf et ' 2 )/2 of the two jets (figure Ct) and of x° bs (figure dl). The 
observable x° bs is defined as {(E — p z )j etl + (E — p z )jet 2 )/ XX — Pz) where tne sum mns over 
all measured particles of the final state. In figures and f the transverse momentum and polar 
angular distributions of the selected tracks are shown. The simulation based on the PYTHIA 
event generator provides a good description of all distributions, after scaling the contributions 
from light quark, charm and beauty events. The scale factors are obtained from fits to the signed 
impact parameter distributions, as described below in section|4| 



4 Quark Flavour Separation 

The fractions of events containing charm and beauty quarks are determined using the same 
method as in previous measurements [5, 6], based on the impact parameter of selected tracks 
which is given by the transverse distance of closest approach to the reconstructed event vertex. 
The signed impact parameter S is defined as positive if the angle between the axis of the asso- 
ciated jet and the line joining the primary vertex to the point of closest approach of the track is 
less than 90°, and is defined as negative otherwise. 

The distribution of S is shown in figure EK- The data are well described by the simulation. 
Due to the long lifetimes of charm and beauty flavoured hadrons the 5 distribution is asym- 
metric, the number of tracks with positive values exceeding the number of tracks with negative 
values. While the component that arises from light quarks is almost symmetric, the c component 
has a moderate asymmetry and the b component shows a marked asymmetry. The asymmetry 
seen at \S\ > 0.1 cm is mainly due to decays of long lived strange particles such as K s . In 
order to reduce the effects of the strange component, tracks with \5\ > 0.1 cm are rejected. The 
significance, defined as the ratio of the impact parameter 5 to its error, is shown in figure |2b f° r 
all selected tracks with \5\ < 0.1 cm. 
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The distributions of Si (figured) and of 5 2 (figure|21i) show the significance of the selected 
track in jets with exactly one selected track associated to the jet (Si) and the significance of 
the track with the second highest absolute significance in jets with two or more selected tracks 
(£2). For jets contributing to the distribution of S2 it is required that the tracks with the first and 
second highest absolute significance in the jet have the same sign of 5. At moderate and large 
values of S 2 the beauty contribution exceeds that from charm. 

In order to substantially reduce the uncertainty due to the resolution of 5 and the light quark 
normalisation the negative bins in the S\ and 5*2 distributions are subtracted from the positive 
ones. These subtracted S\ and S2 distributions are shown in figures and f. The distribu- 
tions are dominated by charm quark events, with an increasing fraction of beauty quark events 
towards larger values of significance. The contribution from light quarks is seen to be small. 

The c, b and light quark fractions in the data are extracted using a simultaneous least 
squares fit of simulated reference distributions for c, b and light quark events, obtained from 
the PYTHIA Monte Carlo simulation, to the measured subtracted S% and S2 distributions (fig- 
ures^ and f). The total number of events before any CST track selection is also used in the fit. 
The Monte Carlo contributions from charm, beauty and light quark events are scaled by factors 
P c , Pb and Pi respectively, which are the free parameters of the fit. Only the statistical errors 
of the data and the Monte Carlo simulation are taken into account. The fit to the complete data 
sample gives scale factors P c = 1.45 ± 0.14, P b = 1.98 ± 0.22, and Pi = 1.44 ± 0.05 and has 
a x 2 /n.d.f. of 13.1/18. The Monte Carlo distributions shown in figures \\\ El El and |7] are scaled 
by these factors. Consistent results are found when alternative methods are used to separate the 
quark flavours, such as the explicit reconstruction of secondary vertices described in section 13 

5 Calculations in Perturbative QCD 

The Monte Carlo simulation programs PYTHIA [18] and CASCADE [20] provide cross section 
predictions in pQCD at leading order. Parton showers are implemented to account for higher 
order effects. PYTHIA uses the DGLAP parton evolution equations [21] while CASCADE 
contains an implementation of the CCFM evolution equations [22] . 

PYTHIA is run in an inclusive mode (MSTP(14)=30 [18]) in which direct and resolved pho- 
ton processes, including heavy quark excitation, are generated using massless matrix elements 
for all quark flavours. The CTEQ5L [23] parton densities are used for the proton and those of 
SaSID [24] for the photon. The charm and beauty quark masses are set to 1.5 and 4.75 GeV re- 
spectively, and the fragmentation is modelled by the Lund string model [25], using the Peterson 
function [26] for the longitudinal fragmentation of beauty and charm quarks. 

Additional Monte Carlo samples of charm, beauty and light quark events are generated using 
the Monte Carlo generator CASCADE with the charm and beauty masses as used in PYTHIA. 
The process •yg — > cc or bb is implemented using off-shell matrix elements convoluted with k t 
unintegrated parton distributions in the proton. In this analysis the parametrisation AO [27] is 
used for the parton distributions. 

QCD calculations to next-to-leading order are performed using the program FMNR [28]. 
FMNR implements the calculation at fixed order in the massive scheme, i.e. charm and beauty 
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quarks are generated dynamically in the hard process via boson-gluon fusion diagrams and the 
parton distributions for the proton and the photon consist only of light quarks (uds) and gluons. 
FMNR provides weighted parton level events with two or three outgoing partons, i.e. the heavy 
quark antiquark pair and possibly a third parton. Values of 1.5 and 4.75 GeV are chosen for 
the c and b quark masses respectively. The renormalisation and factorisation scales are set to 

the transverse masses m t = \Jrn? q + pj M , where pi qq - is the average of the squared transverse 
momenta of the heavy quark and anti-quark. The CTEQ5F3 parameterisation [23] is used for 
the parton distribution functions in the proton. Contributions from processes with resolved 
photons are calculated using the GRV-G HO distributions of partons in the photon [29]. In the 
next-to-leading order prediction these contributions are found to be small (~ 3%). 

In order to compare the parton level calculation with the data, corrections from parton to 
hadron level are applied which are determined using the PYTHIA Monte Carlo event generator. 
The jets at both the parton and the hadron level are reconstructed using the inclusive k t jet- 
algorithm in the p t recombination scheme. The bin-by-bin corrections from parton to hadron 
level are found to be less than ±10% everywhere except in the bins 0.7 <x° bs < 0.85 and 
0.85 <x° bs < 1 where the corrections are about +35% and —15% respectively. 

Theoretical uncertainties of the NLO calculation are estimated by independent variations of 
the renormalisation and factorisation scales by factors of one half and two, and the maximal 
changes to the cross section predictions of 30 — 35% for charm and 20 — 30% for beauty are 
taken as systematic errors. The c (b) masses are varied between 1.4 and 1.6 (4.5 and 5) GeV 
leading to cross section changes of up to ±4%. The cross section variations when using other 
proton structure functions such as CTEQ6M [30], MRSG or MRST1 [31] are less than 8% in all 
regions of the measurement. The latter uncertainty is added in quadrature to the uncertainties 
from the scales and the quark mass. 

6 Cross Section Measurement 

For the measurement of the charm and beauty cross sections, the scale factors P c and P^, which 
are determined from fits of the subtracted significance distributions to the data, are multiplied 
with the cross section predictions of the PYTHIA Monte Carlo simulation. For the measure- 
ment of the differential charm and beauty cross sections the fit is performed separately in each 
bin i. The resulting scale factors P ci and are then multiplied with the bin-averaged cross 
section predictions of the PYTHIA Monte Carlo simulation, divided by the respective bin size. 
In addition, the measured differential cross sections for charm and beauty dijet production are 
divided by the corresponding flavour inclusive cross sections to obtain the fractional contribu- 
tions of events containing charm and beauty quarks. The flavour inclusive dijet cross section 
is measured by correcting the observed number of events before track selections for detector 
efficiencies and acceptances and dividing by the integrated luminosity. 

6.1 Systematic Uncertainties 

Systematic uncertainties of the cross section measurement are evaluated by variations applied 
to the Monte Carlo simulations. The resulting systematic uncertainties of the total charm and 
beauty dijet cross sections are listed in tableland detailed below. 
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Uncertainty [%] 


Source 


Variation 


Charm 


Beauty 


Impact parameter resolution 


©25/xm ©20(Vm 


7 


10 


Jet axis </> direction 


1° shift in (j) 


3 


2 


Track finding efficiency 


2% © 1% 


3 


8 


uds asymmetry 


±50% 


1 


6 


HQ production model (PYTHIA) 


resolved 7, p t dependence 


7 


14 


Fragmentation model 


Peterson / Lund 


1 


2 


Fragmentation fractions 


PDG 


0.5 


1.6 


Hadron lifetimes 


PDG 


0.1 


0.3 


Charged track multiplicities 


MARK-III, LEP, SLD 


1.5 


4 


Jet energy scale 


2% 


6 


5 


Trigger efficiency 




5 


5 


Luminosity measurement 




1.5 


1.5 


Total 




14 


22 



Table 1: Systematic uncertainties of the measured total charm and beauty dijet cross sections. 



• An uncertainty in the 5 resolution of the tracks is estimated by varying the resolution by an 
amount that encompasses the differences between the data and simulation (figureEl). This 
is achieved by applying an additional Gaussian smearing in the Monte Carlo simulation 
of 200 /iin to 5% of randomly selected tracks and of 25 /im to the rest. 

• The uncertainty of the jet axis reconstruction is estimated by shifting the jet axis in az- 
imuth (j) by ±1°. 

• The reconstruction efficiency of central drift chamber tracks is uncertain to the level of 
2% and the efficiency for these tracks to have hits in both r-<p layers of the silicon vertex 
detector is known to 1%. 

• The uncertainty resulting from the shape of the subtracted significance distributions S\ 
and £2 for light quarks (figures and f) is estimated by repeating the fits with the light 
quark Si and S2 distributions varied by ±50% of the default value. 

• The systematic error arising from the uncertainty of the underlying physics model is es- 
timated by varying the contribution from resolved photon processes in the PYTHIA pre- 
diction by ±50%, and by reweighting the p t distribution as predicted by PYTHIA to that 
of CASCADE. These variations lead to changes of the cross sections of ±7% for charm 
and ±14% for beauty. 

• The uncertainties in the description of the heavy quark fragmentation are estimated by 
repeating the fits with Monte Carlo simulation templates in which the Lund Bowler func- 
tion [25] is used for the longitudinal fragmentation instead of the Peterson function. 

• The uncertainties arising from the various D and B hadron lifetimes, fragmentation frac- 
tions and mean charged track multiplicities are estimated by varying the input values of 
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the Monte Carlo simulation by the experimental errors of the corresponding measure- 
ments or world averages. For the fragmentation fractions of c and b quarks to hadrons 
and for the lifetimes of these hadrons the central values and errors on the world averages 
are taken from [32]. For the mean charged track multiplicities the values and uncertain- 
ties for c hadrons are taken from Mark-Ill [33] and for b hadrons from LEP and SLD 
measurements [34]. 

• The uncertainty of the jet energy scale of 2% leads to cross section uncertainties between 
3% at small pf l and 12% at large p{ et and 6% on average, independently of the quark 
flavour. 

• The trigger efficiency is studied using monitoring events from neutral current processes in 
deep inelastic scattering in which the scattered positron triggers the events independently 
of the triggers under study. The uncertainty is determined to be 5%. 

• The luminosity is known to an accuracy of 1.5%. 

Total systematic uncertainties of 14% and 22% are obtained for the measurement of the 
charm and the beauty production cross sections respectively. The total systematic error for the 
flavour inclusive dijet cross section is 8% resulting from the uncertainty of the hadronic energy 
scale (6%), the trigger efficiency uncertainty (5%) and the uncertainty of the luminosity mea- 
surement (1.5%). For the relative contributions of charm and beauty production to the flavour 
inclusive dijet cross section, the statistical errors are added in quadrature and the systematic er- 
rors include those sources that are specific to the charm and beauty cross section measurement. 
The same uncertainties are equally attributed to all bins of the measurement except for the un- 
certainty of the hadronic energy scale for which the uncertainties are determined and applied 
individually in each bin of the measurement. 

6.2 Results 

The total dijet charm photoproduction cross section in the range Q 2 < 1 GeV 2 , 0.15 < y < 0.8, 
p J et K2) > U ( 8 ) GeV and _ Q g < rf* m < 1.3 i s measured to be 

a(ep -> eccX -> ejjX) = 702 ± 67 (stab.) ± 95(syst.)ph. 

For the same kinematic range, the measured beauty cross section is 

a(ep -> ebbX -> ejjX) = 150 ± 17(stat.) ± 33(syst.)pb. 

The predictions from the theoretical calculations are detailed in table |2j The NLO parton 
level calculations are corrected to the hadron level using correction factors as determined from 
PYTHIA (see section^. For charm, the NLO QCD calculation FMNR is somewhat lower than 
the measurement but still in reasonable agreement within the theoretical errors. For beauty, 
FMNR is lower than the data by a factor 1.8, corresponding to 1.6 standard deviations, taking 
both experimental uncertainties and uncertainties in the theory into account. For both, charm 
and beauty, PYTHIA and CASCADE predict a normalisation which is similar to that of FMNR. 
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Charm [pb] 


Beauty [pb] 


Data 
FMNR 
PYTHIA 
CASCADE 


702 ± Q7(stat.) ± 95(syst.) 
5001^ 
484 
438 


150 ± 17(stat.) ± 33(syst.) 

83ig 
76 
80 



Table 2: The measured charm and beauty photoproduction dijet cross sections in the kinematic 

range Q 2 < 1 GeV 2 , 0.15 < y < 0.8, p t Jetl(2) > 11(8) GeV and -0.9 < ^ et ^) < 1.3 in 
comparison to predictions in NLO QCD (FMNR) and from the Monte Carlo programs PYTHIA 
and CASCADE. 



The measured differential cross sections as functions of p x 1 , of 77 and of x° , are shown 
in figure |3] and listed in tableQJ PYTHIA is used to determine the point in the bin at which 
the bin-averaged cross section equals the differential cross section. Both the charm and beauty 
data are reasonably well described in shape. A large difference, however, between the beauty 
data and the NLO QCD calculation is observed in the region of small values of x° bs (figureOf), 
where the prediction lies significantly below the data. In this region, PYTHIA predicts a large 
contribution from events with resolved photons, as indicated in the figures by the dashed-dotted 
lines. PYTHIA describes the shapes of the charm and beauty data distributions, while the 
normalisations are low. The CASCADE prediction is too small in the region of small x° bs , but 
approaches the measurement in the region s° bs > 0.85. 

Differential cross sections are also measured separately for the region x° bs > 0.85 and the 

results as functions of p t jetl and of fj are shown in figureHJand listed in table|U The charm cross 
sections are in good agreement with the NLO QCD calculation both in normalisation and shape. 
The beauty cross sections are also reasonably well described, the agreement being significantly 
better than for the whole range of x° bs . 

The relative contributions from charm and beauty to the inclusive dijet cross sections are 
presented in figure |5] and listed in tables|3] and |4j In figure |5k, the relative contributions are 
shown as a function of x° bs . The data are compared with the PYTHIA Monte Carlo simulation 
which predicts an increase of the relative charm and beauty contributions towards large x° bs 
where direct photon-gluon fusion processes dominate. 

Assuming the charm and beauty quarks to be light, naive quark charge counting predicts a 
value of four for the relative production rates of charm to beauty dijets in direct photon-gluon 
fusion processes. In comparison, the measurement in the region i° bs > 0.85 yields a ratio of 
5.1 ± 1.1 (stat.). In figure|5j5 and c the relative contributions to the dijet cross section are shown 
for the region x° bs > 0.85 as a function of p t jetl and of fj. The ratios are constant within their 
uncertainties. 



7 Heavy Quark Enriched Data Sample 

In order to study the description of the data by the PYTHIA Monte Carlo further, a subsam- 
ple of events is used, in which the fraction of events containing heavy quarks is enriched. For 



11 



this subsample secondary vertices of the heavy hadron decays are reconstructed explicitly using 
events with at least two selected tracks in a given jet. The relative fraction of heavy quarks 
in the event sample is then further enhanced by subsequent cuts on the secondary vertex track 
multiplicity and the decay length significance. For this heavy quark enriched sample the distri- 
butions of the jet transverse momentum, the mean pseudo-rapidity and the observable x° bs are 
presented and the decomposition and shape differences of direct and resolved photon processes, 
as implemented in PYTHIA, are studied. 

Using a method presented in a previous HI publication [5] the primary vertex and the sec- 
ondary vertices of the heavy hadron decays are reconstructed. For each jet the associated tracks 
are used to reconstruct a secondary vertex in an iterative procedure. No definite assignment 
of tracks to vertices is made, but each track is assigned a weight in the range to 1 for each 
vertex candidate, using the weight function of [35]. The smaller the distance of the track to 
the vertex candidate, the larger the weight. For each jet the vertex configuration that minimises 
the global fit x 2 is found iteratively using deterministic annealing [36]. The number of tracks 
which contribute with a weight greater than 0.8 to the secondary vertex is used as a measure 
of the decay vertex track multiplicity. The decay length significance is given by the transverse 
distance between the primary and the secondary vertex divided by its error. The sign is defined 
using the projection of the corresponding vector on the jet direction. Figure |6] shows the de- 
cay length significance distributions for different decay track multiplicities. Good agreement is 
seen between the data and the distributions from the PYTHIA Monte Carlo simulation which 
are scaled by the factors as obtained from the fit to the subtracted significance distributions S% 
and S 2 . Conversely, a simultaneous fit of the decay length significance distributions for the 
different track multiplicities gives results consistent with those from the subtracted significance 
distributions. 

Figure^ shows that in heavy quark events secondary vertices are significantly displaced in 
the direction of the related jet axis. Contributions from light quarks are mainly observed in the 
two-track sample at small decay length significances. In the two-track sample the contribution 
from charm and beauty is similar to that of light quark events while the three and four-track 
samples are dominated by the beauty component, as expected from the mean charged particle 
multiplicity in heavy hadron decays. 

The fraction of events with heavy quarks is enhanced by requiring a secondary decay vertex 
with two or more tracks and a decay length significance larger than 2.0. Applying these cuts, 
the distributions for p t jetl , for fj and for x° bs are shown in figure |7J The data are compared to 
the PYTHIA Monte Carlo simulations. In figures |7^, c and e the decomposition into light, 
charm and beauty quark contributions is indicated. These are determined using the scale factors 
as obtained from the fit to the subtracted significance distributions of the full sample. Good 
agreement between the PYTHIA prediction and the data is seen. The contributions from charm 
and beauty events are 44% and 41% respectively, while 15% of the events contain only light 
quarks. In figures |7J), d and f the same data are shown together with the contributions from 
direct (yg — > qq) and resolved processes, as predicted by the PYTHIA simulation. According 
to the leading order QCD calculation, as implemented in PYTHIA, the fraction of events that 
arise from processes with resolved photons is about 40%. 

In the region x° bs < 0.85 the contribution from resolved processes is enhanced and amounts 
to about 80%. According to PYTHIA the contribution from processes with heavy quark excita- 
tion is by far dominant. In conclusion, the data are well described by the leading order Monte 
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Carlo simulation PYTHIA in which significant contributions from heavy quark excitation pro- 
cesses are predicted. 



8 Conclusions 

Differential charm and beauty dijetphotoproduction cross sections at HERA are measured using 
a technique based on the lifetime of the heavy quark hadrons. The heavy quark cross sections 
are determined using the subtracted impact parameter significance distributions of tracks in 
dijet events. The cross sections are measured as functions of the transverse momentum p t jetl 
of the leading jet, of the mean pseudo-rapidity fj of the two jets and of the observable x° bs . 
Taking into account the theoretical uncertainties, the charm cross sections are consistent both in 
normalisation and shape with a calculation in perturbative QCD to next-to-leading order. The 
beauty cross sections tend to be higher than NLO, by 1.6a for the total cross section, with a 
more significant discrepancy observed in the region of x° bs < 0.85 where processes involving 
resolved photons or higher order contributions are expected to be enhanced. The Monte Carlo 
generator PYTHIA gives a good description of the shape of both the charm and the beauty data. 
The data confirm the expectation within this model that a significant contribution to the heavy 
quark dijet cross section comes from processes with resolved photons. In the region x° bs > 0.85, 
the relative contributions from charm and beauty to the inclusive dijet cross section are found 
to be in agreement within errors with values of 4/1 1 and 1/11, i.e. the naive expectation for the 
direct photon-gluon fusion process, assuming all quarks to be massless. 
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Charm 


p[ 1 range 
[GeV] 


(p? h ) 
[GeV] 


iet 

da/ dpi 1 stat. syst. 
[pb/GeV] 


f cc stat. syst. 


110 14^ 
11. U I'+.J 

14.5 18.0 
18.0 22.5 
22 5 35 


1 9 

1Z. / J 

16.25 
20.0 
27 


1 37 90 90 

49.5 7.6 7.9 

10.6 2.5 1.8 
2 46 86 45 


333 048 044 
V.jjj U.LH-o U.Uff 

0.352 0.054 0.046 
0.270 0.062 0.035 
402 1 40 052 


fj range 


(v) 


da/dfj stat. syst. 
LP D J 


f cc stat. syst. 


-0.90 0.10 
0.10 0.60 

60 1 30 


-0.35 
0.35 

QS 


162 24 22 
674 86 101 
386 S9 61 

JOU Ji Ul 


0.337 0.051 0.045 
0.346 0.044 0.045 

3S8 04Q 047 


x° bs range 




da/dx^ stat. syst. 
[pb] 


f cc stat. syst. 


0.10 0.70 
0.70 0.85 
0.85 1.00 


0.45 
0.77 
0.92 


236 70 34 
1017 271 162 
2994 280 498 


0.250 0.074 0.033 
0.342 0.091 0.047 
0.374 0.035 0.049 






Beauty 


j et i 

p\ range 
[GeV] 


( P r h ) 

[GeV] 


da/dp ! t etl stat. syst. 
[pb/GeV] 


f hb stat. syst. 


11.0 14.5 
14.5 18.0 
18.0 22.5 
22.5 35.0 


12.75 
16.25 
20.0 
27.0 


24.7 6.1 5.5 
9.79 1.87 2.15 
3.37 0.48 0.78 
0.28 0.14 0.07 


0.060 0.015 0.013 
0.070 0.013 0.015 
0.085 0.012 0.018 
0.046 0.022 0.010 


fj range 


(v) 


da/dfj stat. syst. 
[pb] 


/ bb stat. syst. 


-0.90 0.10 
0.10 0.60 
0.60 1.30 


-0.35 
0.35 
0.95 


37.1 5.9 7.8 
112 17 25 

56.2 11.3 12.4 


0.077 0.012 0.016 
0.057 0.009 0.012 
0.052 0.010 0.011 


x" bs range 


« s > 


da/dx^ stat. syst. 
[pb] 


f hb stat. syst. 


0.10 0.70 
0.70 0.85 
0.85 1.00 


0.45 
0.77 
0.92 


44.6 12.8 9.9 
191 60 44 
584 72 135 


0.047 0.014 0.010 
0.064 0.020 0.014 
0.073 0.009 0.015 



Table 3: The measured charm and beauty dijet photoproduction cross sections and the relative 
contributions f cc and / bb to the inclusive dijet photoproduction cross section with statistical and 
systematic errors. 
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Charm 


j et 1 

p[ 1 range 
[GeV] 


(P? h ) 
[GeV] 


da 1 dp ! t etl stat. syst. 
[pb/GeV] 


f cc stat. syst. 


11.0 14.5 
14.5 18.0 
18.0 22.5 
22.5 35.0 


12.75 
16.25 
20.0 
27.0 


73.3 12.7 10.7 
29.0 4.6 4.6 
8.79 1.75 1.46 
1.94 0.70 0.35 


0.348 0.060 0.046 
0.364 0.057 0.048 
0.357 0.071 0.047 
0.455 0.165 0.060 


fj range 




da/dfj stat. syst. 
[pb] 


f cc stat. syst. 


-0.90 0.10 
0.10 0.60 
0.60 1.30 


-0.35 
0.35 
0.95 


115 21 16 
496 70 74 
165 24 26 


0.322 0.058 0.043 
0.462 0.065 0.060 
0.374 0.055 0.049 






Beauty 


j et i 

PI 1 range 
[GeV] 


( P r h ) 

[GeV] 


da/dp ! t etl stat. syst. 
[pb/GeV] 


f hb stat. syst. 


11.0 14.5 
14.5 18.0 
18.0 22.5 
22.5 35.0 


12.75 
16.25 
20.0 
27.0 


15.3 4.3 3.4 
6.19 1.20 1.42 
1.89 0.31 0.44 
0.27 0.18 0.07 


0.073 0.020 0.015 
0.078 0.015 0.016 
0.076 0.013 0.016 
0.063 0.043 0.014 


fj range 


(v) 


da/dfj stat. syst. 
[pb] 


/ bb stat. syst. 


-0.90 0.10 
0.10 0.60 
0.60 1.30 


-0.35 
0.35 
0.95 


29.8 4.7 6.3 
54.4 15.5 12.0 

31.9 5.7 7.1 


0.084 0.013 0.018 
0.051 0.014 0.011 
0.072 0.013 0.016 



Table 4: The measured charm and beauty dijet photoproduction cross sections and the relative 
contributions f cc and f bb to the inclusive photoproduction dijet cross section for the region 
x° bs > 0.85 with statistical and systematic errors. 
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Figure 1 : Distributions of a) the transverse momentum p t 1 of the leading jet, b) the transverse 



momentum p t 2 of the second jet, c) the mean pseudo-rapidity fj of the two jets, d) the ob- 
servable x° bs , e) the transverse momentum of the selected tracks and f) the polar angle of the 
selected tracks. The expectation from the PYTHIA Monte Carlo simulation is included in the 
figure, showing the contributions from the various quark flavours after applying the scale factors 
obtained from the fit to the subtracted significance distributions of the data (see section|U). 
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Figure 2: Distributions of a) the signed impact parameter S of selected tracks, b) the signed 
significance, c) the significance S± of tracks in jets with exactly one selected track, d) the sig- 
nificance S 2 of the track with the second highest significance in jets with two or more selected 
tracks, e) the subtracted signed significance for the sample with exactly one selected track, f) 
the subtracted signed significance for the sample with two or more selected tracks. In b) to f) 
only tracks with an impact parameter \5\ < 0.1 cm are considered. The expectation from the 
PYTHIA Monte Carlo simulation is included in the figure, showing the contributions from the 
various quark flavours after applying the scale factors obtained from the fit to the subtracted 
significance distributions of the data (see section|4|). 
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Figure 3: Differential charm and beauty photoproduction cross sections a,b) da/dp^ 1 , c,d) 
da/dfj and e,f) da/dx^ for the process ep e(cc or bb)X — > ejjX. The inner error bars 
indicate the statistical uncertainty and the outer error bars show the statistical and systematic 
errors added in quadrature. The solid lines indicate the prediction from a NLO QCD calculation, 
corrected for hadronisation effects, and the shaded band shows the estimated uncertainty. The 
absolute predictions from PYTHIA (dashed lines) and CASCADE (dotted lines) are also shown. 
The contribution from resolved processes in PYTHIA (dash-dotted lines) is depicted separately. 
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Figure 4: Differential charm and beauty photoproduction cross sections a,b) da /dpi 1 and c,d) 
da/dfj for the process ep — > e(cc or bb)X — > ejjX in the region x° bs > 0.85. The inner 
error bars indicate the statistical uncertainty and the outer error bars show the statistical and 
systematic errors added in quadrature. The solid lines indicate the prediction from a NLO 
QCD calculation, corrected for hadronisation effects, and the shaded band shows the estimated 
uncertainty. The absolute predictions from PYTHIA (dashed lines) and CASCADE (dotted 
lines) are also shown. 
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Figure 5: Relative contributions from charm (squares) and beauty events (triangles) as a func- 
tion of a) the observable a:° bs , b) the transverse momentum p t jetl of the leading jet for s° bs > 0.85 
and c) the mean pseudo-rapidity fj of the two jets for x° bs > 0.85. The inner error bars indicate 
the statistical uncertainty and the outer error bars show the statistical and systematic error added 
in quadrature. The solid line (shaded area) indicates the absolute prediction from PYTHIA for 
charm (beauty). 
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Figure 6: Decay length significance distributions for samples of events with a secondary vertex 
reconstructed from a) 2, b) 3 and c) 4 tracks. The data (points) are compared with the PYTHIA 
simulation after applying the scale factors as obtained from the fit to the subtracted significance 
distributions of the full sample. 
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Figure 7: Distributions of the heavy quark enriched sample with two or more tracks originating 
from a secondary vertex and a decay length significance larger than 2.0 (see text). The data 
are compared to the PYTHIA simulation for the distributions of a,b) the transverse momentum 
p t jetl of the leading jet, c,d) the mean pseudo-rapidity f} of the two jets and e,f) the observable 
x obs rp^g j £ f t co i umn shows the decomposition of the distribution into beauty, charm and light 
quarks after scaling the PYTHIA predictions by the scale factors obtained from the fit to the 
subtracted significance distributions of the full sample. In the right column the contributions 
from direct and resolved processes in PYTHIA are indicated as shaded histogram and dashed 
line respectively. 
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